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Abstract 

An experimental investigation of a subscale 
HiMAT forebody and inlet was conducted over a 
range of Mach numbers to 1*4. The inlet exhib- 
ited a transitory separation within the diffuser 
but steady state data indicated reattachment at 
the diffuser exit. A finite difference procedure 
for turbulent compressible flow in axisymmctric 
ducts was used to successfully model the HiMAT 
duct flow. The analysis technique was further 
used to estimate the initiation of separation and 
delineate the steady and unsteady flow regimes in 
similar S-shaped ducts 


Introduction 

The inlet designer of advanced tactical air- 
craft must continually expand the stable flight 
and manuaver envelopes. These advanced designs 
often lead to inlets which are characterized by 
short length, sharp curvatures, and large off- 
sets. The associated requirements for high inlet 
performance and inlet/engine compatibility are 
becoming more difficult and close attention must 
be given to the quality of the flow at the dif- 
fuser exit. 

The HiMAT is an example where vehicle design 
considerations have imposed severe geometric con- 
straints on the inlet. In the HiMAT program, the 
test vehicle will demonstrate in flight the inte- 
grated effect of numerous new technologies. 
Although most major components of this RPRV were 
scaled from that of a man-sized high performance 
aircraft system, the inlet's geometric character- 
istics had to be compromised because of aircraft 
volume and engine considerations. The underbelly 
inlet is a fixed geometry, normal shock type with 
a three dimensional cross section. Compared with 
other similar configurations such as the F-16, 
this inlet has significantly shorter length, in- 
creased diffusion and sharper "s" wall curvatures. 

When a subsonic diffuser is shortened, the 
nonuniformity of the flow at diffuser exit in- 
creases end duct separation or unsteady flow can 
exist. Four distinct flow regimes have been de- 
fined for two dimensional diffusers and have been 
related to the geometry of the configuration by 
Kline. 1 These characteristics are shown in 
Fig. 1. Although not defined in the literature, 
similar regimes exist for all other diffusers. 
Since short diffusers of the type used in ad- 
vanced vehicles are prone to separate and produce 
high dynamic distortions at the compressor face, 
the HiMAT diffuser was studied both analytically 
and experimentally, 
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This paper will present results from the 
modeling of the HiMAT inlet using an axisymmctric 
viscous analysis and the subsequent experimental 
subscale HiMAT forebody/inlet distortion test in 
the NASA LeRG 8 by 6 Foot Supersonic Wind Tun- 
nel. The modeling technique was further used to 
generalize the analytical results and delineate 
the steady and nonsteady flow regimes for similar 
S-shaped diffusers. 

Apparatus and Procedure 

Inlet Details 

The RPRV HiMAT was designed for operation to 
Mach numbers in excess of 1*5, The test configu- 
ration was a 0.2884 scale model of the forebody 
and inlet and was tested to Mach 1.4 over a range 
of pitch and yaw angles and a range of engine 
airflows. Figure 2 shows the test model in- 
stalled in the test section of the NASA Lewis 
8 by 6 Foot Supersonic Wind Tunnel. The geome- 
tric shape Of the HiMAT subsonic diffuser is 
shown in Fig. 3 and the axial distribution of 
area in Fig. 4. The relatively constant passage 
height in the initial portion of the diffuser is 
compensated by area changes in the sidewall re- 
gion to provide a smooth ,\rea variation through- 
out tile diffuser* The inlet was coupled to a 
coldpipe-choked plug assembly which maintained 
choked flow conditions a short distance down- 
stream of the compressor face instrumentation 
(Fig. 5). The basic HiMAT inlet does not incor- 
porate vortex generators, but a vortex generator 
study was conducted during the test. The speci- 
fications Of the vortex generators used during 
the study are given in Fig, 6. Vortex generators 
were used on the ramp surface at an x/R 2 of 
5.924 and on the cowl surface at an x/R 2 of 
7.366. The size of the generators was based on 
computed flow conditions on the ramp surface at 
the upstream generator location. No adjustments 
in generator size were made for the side wall or 
for the downstream cowl surface flow conditions. 


Instrumentation 

The primary total pressure instrumentation 
is shown in the cutaway sketch (Fig. 5). Four 
boundary layer rakes were placed in the diffuser 
m the ramp surface to detect the predicted flow 
tfcparation. The upstream rake was placed ahead 
the anticipated separation whereas the other 
chcee rakes were used to define the separation 
region. The rakes were staggered circumferen- 
tially to avoid interference effects on down- 
stream instrumentation. Forty combination steady 
state-dynamic total pressure probes were used at 
the compressor face to measure the diffuser exit 
conditions. These probes were arranged in a 5 
ring by 8 rake probe configuration. The pressure 
probes were located at the center of equal 
areas. Only the fluctuating component of pres- 
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eure van recorded from each dynamic probe. The 
fluctuating component of each pressure transducer 
was recorded on FM multiplexed tape. This compo- 
nent was later added to the corresponding time 
averaged (steady state) pressure during data 
analysis. The frequency response of the dynamic 
compressor face probes was flat to about 
2000 Hr.. Ail analysis of the dynamic data was 
performed on DYtfADEC by AFFDL. The local static 
pressures were measured on both the ramp and cowl 
centerlines and by a row of side wall static tops. 


Results and Discussion 

The short length, high diffusion and sharp 
8-curvature of the HiMAT subsonic diffuser sug- 
gested a possibility of flow separation. 7. 
evaluate this possibility, the finite-difference 
procedure of Ref. 2, for turbulent compress*' hie 
flow in axisymmetric ducts was used. This analy- 
sis procedure first obtains a potential flow 
solution which is used to establish a generalized 
orthogonal curvilinear coordinate network. The 
turbulent viscous effects are then retained in 
the formulation of the problem and the analysis 
is reinitiated to obtain a time marching solution 
with viscous effects Included throughout. The 
program approximates the flow blockage and fric- 
tional effects of struts in the analysis. 

To model the conformal inlet by an axi- 
symmetric code, the geometry was simulated by a 
segment of an annular diffuser using the proce- 
dure of Ref. 3. The geometric simulation was 
constrained to preserve: (1) the HiMAT diffusor 

area ratio, (2) the HiMAT entrance height at the 
centerline, (3) the HiMAT exit height at the 
centerline (the compressor face diameter), and 
(4) the HiMAT offset on the ramp Surface. These 
conditions also ensure that the offset on the 
cowl surface is maintained. The four constraints 
determine the exit and entrance radii for the 
inner and outer contours of the equivalent annu- 
lar geometry. Struts were added of an appropri- 
ate number and thickness to provide the proper 
entrance area of HiMAT, The ramp and cowl con- 
tours of HiMAT at the centerline are then also 
preserved. The resultant geometric simulation is 
schematically shown in Fig. 7 and the comparison 
of the modeled and actual diffuser cross sections 
are compared in Fig. 8. This simulation does not 
accurately model the areas except at the diffuser 
entrance and exit and does not model secondary 
flows within the diffuser. The vehicle operating 
conditions chosen for the analysis were a free 
stream Mach number of 0.9, an altitude of 30 000 
ft, a throat Mach number of 0,7 and zero angle- 
of-attack and yaw. An assumed boundary layer 
flow blockage of 0.9% was used. Figure 9 shows 
the resulting velocity ratio plot and Fig. 10 
shows the resulting wall shear stress distribu- 
tions. Kline in Ref. 1 has suggested that high 
subsonic inlet Mach numbers tend to promote early 
transitory unsteady flows in short conical dif- 
fusers. Since the entrance Mach number of this 
diffuser is also high subsonic and the analysis 
does not include secondary effects, the shear 
distribution on the ramp surfac? suggests a high 
probability that separation will occur at an 
x/Rrji of about 0.55. No collusion can be 
drawn about the effects of this separation on the 
peak instantaneous distortion at the compressor 
face or whether reattachment will occur within 
the diffuser. The scale model of the HiMAT in- 


let, including the forebody and canards, was 
therefore built and tested in the tewis 8x6 SWT 
to determine the separation characteristics and 
the associated effects on pressure recovery and 
distortion levels. 

The axial distributions of static pressure 
on the ramp, cowl, and sidewall arc given in 
Fig. U* The distributions shown are for zero 
angle-of-attack and yaw and a free stream Mach 
numHr of 0.9. The pressure distributions show a 
strong expansion on the cowl and side surfaces 
near the diffuser entrance. A small expansion 
occurs on the ramp walls downstream to a x/R c 
of about 8.4. This expansion is due to local 
boundary layer growth and wall curvature ef- 
fects. The results of the analysis are also 
shown in the figure. The separation appears to 
take place near the sudden change in slope of the 
wall static pressure at X/R 2 p 10.3 shown in 
Fig. 11. 

The flow along the ramp surface was surveyed 
at four axial locations. The total pressure dis- 
tributions ore shown in Fig. 12. The wall static 
pressures as obtained by interpolation from the 
ramp distribution of Fig. 11 are included on the 
profiles. The profile from the upstream rake is 
typical of a turbulent boundary layer and shows 
little free stream recovery loss. Each of the 
downstream rake profiles indicate high recovery 
loss and r low separation and an attendant loss in 
recovery. The initiation of the separation is 
therefore between an axial location of x/R c « 

9, 503 and 10.939. 

The steady state total pressure profiles at 
the compressor face are shown in Fig. 13, The 
recovery is plotted for each of the four quad- 
rants. The local wall static pressure ratios are 
also shown. These data were measured with the mid 
diffuser rakes removed to eliminate any possible 
rake wake effect. Recovery on the cowl surface 
was generally very good except in the near wall 
region# The recovery profiles on the ramp sur- 
face clearly show a large reduction in pressure 
on the top centerline of the diffuser which was 
downstream of the flow separation. The near wall 
total pressures are higher than the measured wall 
static pressures. This indicates that flow re- 
attachment has occurred. Examination of the dy- 
namic pressure however indicates that transitory 
separation is present. A typical time history of 
the dynamic component of total pressure obtained 
from a typical dynamic pressure probe at the com- 
pressor face is shown in Fig, 14. The pressure 
trace shows a one-sided character. When the mean 
total pressure 1 b assigned the steady state level 
of pressure as independently measured, the down- 
ward Bpiking extends from this mean steady state 
total pressure level down to the steady state 
static pressure level. This strongly suggests 
that at those instants of time when downward 
spiking occurs in the pressure traces, a transi- 
tory separation occurs at the compressor face. 
Additional creditability to this argument is ob- 
tained by noting that the upper level of the 
envelope of time varying pressures corresponds 
approximately to the mean free stream total pres- 
sure. The downward spiking in the pressure sig- 
nals occurred randomly and was not associated 
with any particular frequency. The spiking were 
usually noted on the probes of the outer ring and 
on the probes downstream of the diffuser ramp 
separation. On occasion, however, they were 
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noted on each of the other dynamic pressure 
probes. 

One well known technique for eliminating 
duct flow separation is the use of vortex genera- 
tors* The vortex generators described previously 
were installed upstream of the separation zone to 
energize the low energy fluid near the wall. The 
effectiveness of the vortex generators was deter- 
mined by using the downstream diffuser boundary 
layer rake* The total pressure profile is com- 
pared to the comparable profile without genera- 
tors in Fig. 15* The profile shows that separa- 
tion was eliminated by the vortex generators. 
Examination of the dynamic traces confirmed that 
the transitory separation also had been elimi- 
nated. 

The average total pressure recovery measured 
at the compressor face is shown in Fig. 16. Re- 
sults are presented for operation both with and 
without vortex generators, The data shown is for 
nominal vehicle operating conditions which ares 
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The recovery which corresponds to the zero 
angle-of-attack data presented in Figs* 11 to 15 
is shown by the filled symbol. The re cove rys are 
very good over the entire range of ope vat ion de- 
spite the presence of the upstream separation on 
the diffuser surface and the compressor face 
transitory separation* The addition of the vor- 
tex generators reduced the recovery about 1 per- 
cent over the Mach number range. 

The distortion in terms of APRS n the dis- 
tortion methodology for the J85-GE-21, is pre- 
sented in Fig, 17. The distortion is the peak 
instantaneous distortion obtained from the analy- 
sis by AFFUl using DYNADEC. Each value of dis- 
tortion was obtained from a 30 second recording 
of dynamic pressure data, The data was processed 
using a 500 Hz low pass filter to account for 
engine response. The recommended limit for the 
J85-GE-21 engine is also shown for operation at a 
40 000 foot altitude. This limit assumes a nor- 
mal engine operating with a unrestricted throt- 
tle. The distortion delivered by the inlet is 
lower than the engine limits for all operating 
conditions shown. The figure also shows the dis- 
tortion when vortex generators were added. The 
peak distortion is significantly reduced when the 
vortex generators were installed except at the 
low Mach numbers. The generally higher distor- 
tion without the vortex generators is a result of 
a difference in the character of the distortion, 
The basic pressure patterns and hence the distor- 
tion changes when vortex generators are installed 
in the diffuser. Without generators, the peak 
distortion always occurs during that instant of 
time when transitory separation is present at the 
compressor face. The pressure distribution is 
characterized by a large "hole" downstream of the 
diffuser separation. This gives rise to large 
values for both the circumferential and radial 
distortion indices. The J85-GE-21 engine can 
however tolerate higher levels of circumferential 
distortion when accompanied by elevated levels of 
radial distortion. The resultant measured in- 


stantaneous distortion in terms of & 3 RS is mod- 
est* When generators are installed the large 
"hole 11 in pressure is no longer present* As a 
result the radial and the circumferential distor- 
tions are reduced to about half the previous val- 
ues. The value of <6PRS is not reduced by the 
same amount however since the reduction in radial 
distortion causes an increase in £PRS which is 
not offset by the reduced circumferential distor- 
tion. 

The experimentally measured profiles of 
total pressure at the diffuser exit are compared 
to the predicted profiles in Fig* IB. In the 
modeling no engine hub was included* The experi- 
mental profiles that ore presented are for the 
centerline rakes on the ramp (top) and cowl (bot- 
tom)* The ramp side measured pressures do not 
agree with the predictions because of the effects 
of the separation. The agreement between the 
predicted and experimental pressure profiles on 
the cowl side of the diffuser is considered ex- 
cellent. 

The agreement shown between the static pres- 
sures (Fig. II) and total pressures (Fig. IB) 
suggest chat the modeling technique may be gener- 
alised to define or delineate the separation 
characteristics of similarly shaped diffusers. 
Kline in Ref. 1 showed that the flow regimes of 
conical diffusers and two dimensional diffusers 
correlate in terms of 0, the diffuser half angle 
and the nondimensional length. Specifically, the 
area ratio can be expressed as 

A 2 /Ai « 1 + 2 (sin £)L/L ¥ef 

For two dimensional diffusers the reference 
length used in the correlation was determined to 
be the inlet passage height and for conical dif- 
fusers it was the inlet radius. The results of 
Kline have been recast in the form presented in 
Fig. 19. The "line of first stall 1 ' and the "line 
of appreciable stall" as determined by Kline are 
shown. The area ratio presented is a measure of 
the diffusion whereas the mean wall angle is a 
measure of the flow turning. Kline’s correlation 
in this form suggests that it can be used as Che 
basis of defining the flow for HiMAT type dif- 
fusers. 

The modeling procedure previously presented 
was used to analyize the HiMAT type diffusers, A 
matrix of configurations was generated by varying 
the area ratio (constraint no* 1) and the axial 
length while preserving the other HiMAT Geometric 
constraints. Struts for the axisymmetric dif- 
fuser segment analysis were selected to preserve 
the diffuser exit area. This matrix of geome- 
tries is considered to be HiMAT-simiiar- 
geometries since the ramp and cowl centerline 
contours are preserved. The axial length was 
varied by applying a scale factor to the axial 
coordinates. For a given area ratio, configura- 
tions with various lengths were then analyized 
until a diffuser separation was analytically en- 
countered. The results of the analysis are pre- 
sented in Fig, 19 in terms of the area ratio and 
the average wall angle of the surface on which 
separation was obtained with the wall shear 
stress as a parameter. The value of wall shear 
stress presented is the minimum value within the 
diffuser. The line of zero wall shear stress 
agrees very well with the line of appreciable 
stall. The second set of conditions are for a 
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wall shear coefficient of Q»0002» This value was 
obtained in the HiMAT analytical results of 
Fig. 10 and was shown experimentally to result in 
transitory separation. These results are in good 
agreement with Kline's line of first stall* 

The agreement between the HiMAT class of 
diffusers and Kline's Stanford two dimensional 
diffuser correlation is considered adequate for 
defining flow regimes. 

It should finally be noted that all of the 
results presented in this paper are for HiMAT 
rnmp and cowl centerline contours end the appli- 
cability of the results to other configuration 
has not been determined. 


Concluding Remarks 

A wind tunnel investigation of the HiMAT 
forebody and inlet was conducted over a range of 
Mach numbers to 1.4. The experimental results 
were compared to analytical predictions from a 
finite difference procedure for turbulent com* 
pressible flow in axisymmetric ducts. The re* 
suits indicate the fol lowing i 

1. The average total pressure recovery at 
the compressor face was high over the nominal 
vehicle operating conditions. The peak inetanta* 
neous distortion obtained during operation at 
nominal conditions was less than the recommended 
limit for the J85-GE-21 engine. 

2* Separation was noted for all operating 
conditions within the subsonic diffuser. Steady 
state total pressure profiles at the diffuser 
exit indicated reattached flow; however, dynamic 
pressure data indicated the presence of a transi* 
tory separation. The addition of vortex genera- 
tors in the diffuser eliminated the duct separa- 
tion. 

3. The duct flow was analytically modeled as 
a segment of an axisymmetric annular diffuser. 

The shear stress results indicated a high risk of 
duct separation. The static pressure predictions 
on the walls of the diffuser and the total pres- 
sure distributions at the diffuser exit compared 
favorably with the experimental da* in those 
portions of the flow field not af*,cted by the 
flow separation* 

4. The analysis technique was used to esti- 
mate the initiation of separation and delineate 
the steady and unsteady flow regimes in similar 
s-shaped ducts. 


Symbols 


full scale tip radius at inlet throat 
for annular diffuser analytical 
program (214.80 cm) 
full scale compressor face radius 
(22. 987 cm) 

x axial distance 

y vertical distance 

Subscripts; 

I local 

ref reference 

T tip or cowl for annular diffuser analyt- 

ical program 

II hub or ramp for annular diffuser analyt- 

ical program 

0 free stream 

1 diffuser entrance station 

2 diffuser exit station 
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A area 

d distance from surface 

h vertical height 

H annular height at diffuBer exit station 

(tip radius-hub radius) 

L length 

L ref reference length 

p static pressure 

P total pressure 

r radial distance from diffuser exit 

centerline 
R radius 
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figure 1. - Mom regimes for two-dimensional 
diffusers iref. II. 


Figure 2 . - Model installation in the N^A Lewis 8- by 6-foot Supersonic Wind Tunnel 
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Figure 3. - Inlet Internal contours, 
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Figure 4. - Diffuser area distribution. 






Figure 7. - Sketch of annular diffuser segment used In the analysis 
program to simulate the HIMAT diffuser. 
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Figure 10, - Wall shear stress distribution, 
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Figure 11. - Analytical and experimental distribu- 
tions of wall static pressure, 
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Figure 13. - Local distribution of compressor face recovery. 
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Figure 14, * Time history of typical compressor face total pressure. 
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Figure 15. - Subsonic diffuser 
total pressure profile with and 
without vortex generators In* 
stalled at diffuser station, 
x/R 2 , 12.265, 
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Figure 16. - Compressor face pressure recovery with 
and without vortex generators. 
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Figure 17, - Compressor lace Instantaneous distortion 
with and without vortex generators. 



Figure 18. • Comparison ol analytical 
and experimental compressor (ace 
total pressure profiles. 
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Figure 19. • Dilfusor flow regimes, 
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